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In this paper we review basic properties of coupled periodic and aperiodic nanopillar waveguides. A coupled 
nanopillar waveguide consists of several rows of periodically or aperiodically placed dielectric rods (pillars). 
In such a waveguide, light confinement is due to the total internal reflection, while guided modes dispersion 
is strongly affected by the waveguide structure. We present a systematic analysis of the optical properties 
of coupled nanopillar waveguides and discuss their possible applications for integrated optics. 
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1 Introduction 

Photonic crystals (PhCs) are known for offering unique opportunities for controling the flow of light by 
acting as waveguides, cavities, dispersive elements, etc [1-4]. Photonic crystal waveguides (PCW) are 
one of the promising examples of PhCs applications at micron and sub-micron length-scales. They can be 
formed by removing one or several lines of scatterers from the PhC lattice (Fig.[T^). PCW based on PhCs 
with different two-dimensional (2D) lattices of both air holes in a dielectric background and dielectric rods 
in air were reported [1-3]. Light confinement in PCW is obtained due to a complete photonic bandgap 
(PBG), in contrast to the standard guiding mechanism in a conventional dielectric waveguide (Fig.[T}3). It 
was theoretically predicted that a PhC waveguide can possess loss-free propagation as soon as a guiding 
mode falls into a complete PBG. However, progress in PhC research has revealed that losses are inevitable 
and sometimes might be rather high even in spite of broad PBG. Special optimization efforts are now 
intensively applied for decreasing optical losses and the results are quite promising [5,6]. 

At the same time, PBG guiding is not the only waveguiding mechanism in a PhC. Unique anisotropy 
of PhCs can cancel out the natural diffraction of the light, leading to the self-guiding of a beam in a non- 
channel PCW [7-9]. The common principle of index guiding (guiding due to total internal reflection) can 
be also found in periodic systems. It is rather straightforward if a waveguide is organized as a defect in 
a lattice of holes in a dielectric material. Then, the channel itself has higher index of refraction than the 
average index of the drilled or etched medium. Topologically inverted systems like periodic arrays of rods 
or nanopillars placed in air can also provide waveguiding due to index difference [2]. However, fabrication 
of rod arrays on the nanoscale is a relatively difficult technological problem. 

The recent progress in the fabrication of alternative structures, namely nanopillar waveguides (NPW), 
has proved the relevance of their study not only as a useful theoretical model. For example, two-dimensional 
(2D) silicon-on-insulator (SOI) pillar PhC have recently been fabricated and characterized [10]. Sandwich- 
like structures have also been successfully realized in GaAs/AlxOy material system [11]. Membrane-like 
structures have been realized, based on polymer membranes incorporating Si rods [12]. Recently, various 



Con-esponding author: e-mail: chigrin@th.physik.uni-bonn.de, Phone: +49 228 73 2046, Fax: +49 228 73 3223 



pss data will be provided by the publisher 



2 



D. N. Chigrin, S. V. Zhukovsky, A. V. Lavrinenko, and J. Kroha: Coupled nanopillar waveguides 




Fig. 1 Optical waveguides: (a) photonic crystal waveguide, (b) dielectric waveguide, (c) nanopillar waveguide, and 
(d) coupled nanopillar waveguide. 

combinations of active materials inserted in single nanowires or arrays of nanopillars have been under at- 
tention as well [13]. It is important to point out that all of the above mentioned studies do not only present 
a successful practical realization of the pillar PhC structures, but also report transmission efficiencies and 
out-of -plane radiation losses comparable with the 2D PhC based on hole geometry. 

A one-dimensional (ID) chain of rods placed at equal distance from one another (Fig.[T];) possesses 
guiding properties as was shown by Fan et al. [2]. The fundamental mode of such a periodic nanopillar 
waveguide lies below the light line and below the first PBG corresponding to the 2D PhCs with a square 
lattice of the same rods. Guiding is due to total internal reflection. A better confinement of light can be 
achieved, if several ID periodic chains are placed in parallel (Fig.[TJl) [14]. Such waveguides are called 
coupled nanopillar waveguides (CNPWs) and are designated as Wn, where n is the number of parallel 
rows comprising the CNPW. In building a CNPW both the longitudinal and the transverse relative shift 
between individual waveguides can be arbitrary, and thus, a high flexibiUty in dispersion engineering can 
be achieved. 

In this paper, we review basic properties of coupled nanopillar waveguides and discuss their possible 
applications for integrated optics. In Section[2l a CNPW is introduced and possible ways to tune the CNPW 
dispersion are discussed. The transmission efficiency of 2D and 3D CNPWs is reported in Section [3] The 
route to improve the coupling between a nanopillar waveguide and an external dielectric waveguide (like 
an optical fiber) is discussed in Section|4]with respect to aperiodic NPWs. Possible applications of coupled 
periodic and aperiodic nanopillar waveguides are discussed in Section|5] Section|6]concludes the paper 



2 Dispersion engineering 

2.1 Dispersion tuning 

In reference [2] it was shown that a single row of periodically placed dielectric rods is effectively a single- 
mode waveguide within a wide frequency range (Fig. |2] left panel). It has a well confined fundamental 
mode. Attaching one, two or more identical Wl waveguides in parallel to the original one produces a 
coupled-waveguide structure [14]. It is well known in optoelectronics that this leads to the splitting of the 
original mode into n modes, where n is the number of coupled waveguides [15]. 

In Fig. |2] dispersion diagrams for Wl, W2, W3 and W4 CNPWs are shown. All rods are placed at 
the vertices of a square lattice. To model a CNPW dispersion we used the plane-wave expansion method 
(PWM) [16]. The supercell consists of one period in the z direction and 20 periods in the x direction, where 
n periods occupied by dielectric rods were placed in the center of the supercell. The waveguide is oriented 
along the z-axis (Fig.|2]l. The calculations were performed for 2D structures and for TM polarization. The 
n modes of the CNPW are bound between the T — X and X — M projected bands of the corresponding 
infinite PhC of a 2D square lattice of rods (Fig.|2] dashed lines) [14]. All modes are effectively localized 
within the waveguide region. Near the irreducible Brillouin zone (IBZ) boundary the dispersion is strongly 
affected by the system periodicity. 
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Fig. 2 Dispersion diagrams for CNPWs with 1, 2, 3 and 4 rows. The insets show a sketch of the waveguides. The 
coordinate system together with the first quadrant of the first Brillouin zone of the square lattice with period a are also 
shown. The grey areas depict teh continuum of radiated modes lying above the light line. Guided modes are shown as 
black solid lines. The projected band structure of the infinite 2D PhC is shown as dashed lines. Here e = 13.0 and r = 
0.15a. 



It is well known that by varying the filling factor, i.e. the rod radius, and the dielectric constant of the 
rods, one can tailor the frequency range and slope of the PhC bands. Taking into account that the CNPW 
modes are bound hy T — X and X — AI bands of the corresponding infinite 2D PhC, a proper frequency 
adjustment of nanopillar waveguide modes can be done by changing these two parameters. Decreasing the 
dielectric constant of the rods, while keeping their radius constant, pushes the bundle of n CNPW modes 
to higher frequencies. The modes shift towards lower frequencies, if the nanopillar radius increases, with 
fixed dielectric constant. In general, the mode tuning follows the rule: the larger the average refractive 
index of the system, the lower the mode frequencies [14]. 

Another option for tuning the mode dispersion of CNPW is to change the distance between individual 
waveguides, the transverse offset. Examples are shown in figure [3] for two transverse offsets, d = 0.5a 
(left) and d — 2.0a (right). In these cases the rods are situated at the vertices of a rectangular lattice. 
While the mode overlap of individual waveguides is larger (smaller) for close (far) positioned waveguides, 
the coupling strength is stronger (weaker). For two identical waveguides, this in turn results in stronger 
(weaker) mode splitting, (3± = (3 ± k, with respect to the propagation constant (3 of the uncoupled NPW. 
Here k is a coupling coefficient [15]. Note, that the CNPW mode frequencies are still bounded by the 
position of the projected band structure of the corresponding infinite rectangular PhC (Fig. [3] dashed Unes). 



0.5 u H h > ' 0.5 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



Wave vector (2jr/a) 

Fig. 3 The same as in Fig.|2]for CNPWs with different transverse offsets, d — 0.5a (left) and d = 2.0a (right). The 
insets show a sketch of waveguides and coordinate system. 
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The last parameter which may affect the dispersion of a CNPW is the longitudinal shift between its 
individual rows. In figure |4] the dispersion diagrams for CNPW with rods placed in the vertices of a 
triangular lattice are shown for W2, W3, W4 and W5 waveguides. The orientation of the waveguides 
coincides with the T — X direction of the triangular lattice. The mode splitting in a "triangular lattice" W2 
waveguide strongly depends on the propagation constant (Fig. IH left panel), being large for small /3 and 
vanishing near the IBZ boundary. This is in contrast to a "square lattice" W2 waveguide (Fig. |2]i, where 
the mode splitting is approximately constant for all propagation constants. The mode multiplicity near 
the IBZ boundary leads to regions with negative dispersion (backward propagating waves) of the second 
mode. For CNPWs with the number of rods larger than two (Fig.|4] right panels) it results, furthermore, in 
the formation of mini-bandgaps and multiple backward waves regions in the dispersion. Note that in spite 
of the complex nature of the mode splitting, CNPW modes are still bounded by the projected bands of the 
corresponding triangular lattice PhC. 

The longitudinal shift 6 can be arbitrarily set to any value between 6 — and S = 0.5a. The concomitant 
dramatic changes in the CNPW dispersion are illustrated in figure|5]for the case of W2 waveguide. Starting 
from the simple mode splitting for 5 = one can have a very flat second band for S w 0.25a, with negative 
dispersion regions in the second band for 5 > 0.25a and degenerate first and second bands at the IBZ 
boundary for 6 — 0.5a. By combining such shifted W2 waveguides and appropriately chosing the rod 
radius and transverse offset one obtains large flexibility in designing CNPWs with anomalous dispersion 
in the frequency range of interest. 



2.2 Coupled mode model 

To understand qualitatively the physical mechanism of the anomalous dispersion presented in the last 
example (Fig.jS), the coupled mode theory (CMT) can be used [17]. Being an approximate theory, CMT 
nevertheless manages to combine a simple physical model with accurate qualitative and even quantitative 
results [18]. In what follows, two identical coupled periodic waveguides a-Wl and b-Wl are arranged in 
a W2 CNPW. The second waveguide, b-Wl, is shifted by 6 with respect to the first one (Fig.|5] inset). We 
limit ourselves to the scalar CMT, which in our case corresponds to the TM polarization. 

The modes of the W2 waveguide are defined as the solutions of the 2D scalar wave equation 

/ 52 ^2 \ 

( + ^ ) ^ (2^' ^) + i^, ^) E (x, z) = 0, (1) 

where the dielectric function of the composite structure is simply the sum of dielectric functions of the 
two Wl waveguides, e (x, z) = £a [x, z) + Sb {x, z). Here ko — uj / c is a wave number in vacuum. We 
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are looking for a solution of Eq. ^ in the form of a linear combination of the propagating modes of two 
isolated Wl waveguides [19], which allows us to separate spatial variables in the form 



E (x, z) = (x) {fa (z) e^'^' + ba (z) e'^') + {x) e-'P' (/, (z) e'^^^ + h [z) e'^') . (2) 

Here (z) = Fm {z) exp {i {(3 — (3q) z) and bm {z) — B,n {z) exp (— i (/? — z) are the slowly varying 
amplitudes of forward and backward propagating modes near the Bragg resonance condition of a single 
periodic Wl waveguide with period a and /3o — Ti/a. The functions 'i'a {x) and iP'b [x) represent the 
transverse field distributions, and indexes m — a,b refer to a-Wl and b-Wl waveguides, respectively. The 
spatial shift between the two Wl waveguides is accounted for by the corresponding phase shift e'"^^^ of 
the field of the b-Wl waveguide. Here (3 is the propagation constant of a homogenized Wl waveguide of 
width / = 2r and dielectric constant £cfr (x) — ° dze {x, z). The dependence of the propagation 

constant (3 on frequency is given by the standard planar waveguide dispersion relation [15] 




where we have introduced the effective index of refraction of the homogenized waveguide n^a = ^/e^. 
The transverse field distributions Wa [x) and ^i, [x) obey the scalar wave equations 

- ^'m {x) + klsQm {x) tf,„ {x) ^ 

with m — a,b and the transverse dependent dielectric functions £om (x) being a z-average dielectric con- 
stant of the m-th waveguide. Substituting the mode expansion (|2]) into the scalar wave equation ^ and 
expanding the dielectric constant e {x, z) in a Fourier series with respect to z, 

e {x, z) = e„o {x) + {x) + ^ {sai {x) + su {x)) g-^'C^V^.)^^ (4) 

with the Fourier coefficients Smi {x), one obtains after some lengthy but straightforward derivations a sys- 
tem of four ordinary differential equations relating slowly varying amplitudes of the forward and backward 
propagating modes in the two Wl waveguides. 
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Wave vector (27r/a) 

Fig. 5 Dispersion diagrams of the W2 CNPW for different longitudinal shifts 5. The insets show a sketch of two 
parallel periodic waveguides with transverse offset d and longitudinal shift 5. Here e — 13.0, r = 0.15o and d = a. 
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Wave vector (27r/a) Wave vector (27r /a) Wave vector (27r/a) 

Fig. 6 Dispersion diagrams of two coupled periodic waveguides (solid lines) within the framework of coupled mode 
theory for three values of the longitudinal shift 5 = 0.0 (left), S = 0.5 (center) and 5 — 0.25 (right). The dotted 
line is the dispersion of a homogenized waveguide folded into the first Brillouin zone. The dashed lines are the folded 
dispersions of two coupled homogenized waveguides. Here rieff = 1.5, 1 = 0.3a, kq = 0.06 and Ka ~ 0.03. 



For the propagation constant close to the Bragg point /Sq ~ ir/a, the system matrix AI has the form 



M = 



Vo 







-Kq 







P- 



Po 



(6) 



To simplify the following analysis, we have kept only two coupling constants, namely kq, accounting for 
the coupling between two homogenized waveguides, and Kq, describing the waveguide's intrinsic periodic 
structure. These coupling constants are defined in a usual way, as overlap integrals of the transverse field 
distributions with the corresponding Fourier coefficients of the dielectric function expansion. The resulting 
propagation constants of the supermodes of the W2 waveguide are given as the eigenvalues of the system 
matrix AI ^ 

/3w2H =/3o + A/3H, (7) 



with 



A/3 (w) = ±y iP (uj) - Pof + «g - ^2 ± «oy^4 (/? ~ /3o)' - 2^1 + 2kI cos {2ttS). (8) 

The implicit dependence of the propagation constant on frequency is given via the dispersion relation (|3]l 
of a planar homogenized waveguide. 

In figure |6]the dispersion diagram of W2 waveguide calculated using Eqs. (|7]|8]l is presented for three 
values of the longitudinal shift 6 — 0.0 (left), 6 — 0.5 (center) and S — 0.25 (right). See figure caption for 
further details on the parameters. The dotted line shows the dispersion of a planar homogenized waveguide 
calculated using the dispersion relation (O and folded back into the first Brillouin zone by the Bragg wave 
vector corresponding to the periodic Wl waveguide. By setting the self-action coupling constant, Ka, to 
zero and choosing some finite value for the inter-row coupling constant, kq, one can reproduce the simple 
band splitting within the CMT model. The split modes are shown as dashed lines (Fig.|6]l. 

To analyze the influence of the periodic structure and the longitudinal shift on the split band structure, 
we first consider zero longitudinal shift, S ~ 0.0. In this situation the detuning of the propagation con- 
stant from the Bragg wave vector, Pq, is given by A/3 = ±^J (A ± kq)^ — where A = {P — Po) 
is the detuning of the propagation constant of the homogenized waveguide from the Bragg point. The 
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Fig. 7 Left: Dispersion diagram and transmission spectra of a 2D "square-lattice" W4 CNPW. The fundamental 
mode was excited. Here e = 13.0 and r — 0.15a. Right: Dispersion diagram and transmission spectra of a 2D W3 
"triangular-lattice" CNPW. Solid line - even excitation, dashed line - odd excitation. Here e — 13.0 and r = 0.26a. 



propagation factor of the supermodes is given by the exponential e^^'^o^e which cor- 

responds to propagating modes only if (A ± kq)'^ — > 1. In the opposite situation, there are two 
bandgaps at Bragg wave vector /?o with central frequencies corresponding to A = zLkq. These bandgaps 
are due to the destructive interference of the first forward propagating and the first backward propagating 
supermodes and the second forward propagating and the second backward propagating supermodes, re- 
spectively, as can be seen from the left panel of figure|6] In the case of half-period shifted Wl waveguides, 
S = 0.5, the detuning of the propagation constant and the supermode propagation constant are given by 

A(3 = ± ^-^/A^ — kI± kq^ and e.±^(Po±f^o)z g^^{V^~^) ^ jjj ^^lis case two bandgaps exist at shifted 
Bragg wave vectors j3o ± with central frequency at A = 0.0. This corresponds to the destructive 
interference of the first forward propagating and the second backward propagating supermodes and vice 
versa (Fig.|6] center). It is important to mention here that at the Bragg condition, /3o, (IBZ boundary) the 
first forward propagating and the first backward propagating supermodes are in phase, which leads to the 
degeneracy of the first and second bands at the IBZ boundary (Fig. |4]l. The shift of the Bragg condition 
away from the IBZ boundary, /3o ± kq, is a reason for the appearance of a region with negative dispersion 
in the second band (Fig.|4|i. In the case of arbitrary shift between the Wl waveguides a destructive interfer- 
ence takes place between all possible pair combinations of forward and backward propagating supermodes 
leading to the formation of four bandgaps and anomalous dispersion (Fig. |5]l. In the right panel of fig- 
ure |6] an example of a CMT dispersion diagram is shown for the case of quarter-period shifted waveguides, 
S = 0.25. 



3 Transmission efficiency 

An important characteristic of the novel waveguides is their transmission efficiency. To analyze transmis- 
sion efficiencies of different 2D and 3D CNPW, the finite difference time domain (FDTD) method [20] 
with perfectly matched layers as absorbing boundary conditions at all sides and a resolution of 16 grid 
points per lattice constant is used here. The modes are excited by a Gaussian-shaped temporal impulse, the 
Fourier transform of which is broad enough to cover the frequency range of interest. Fields are monitored 
by input and output detectors. The transmitted wave intensities are normaUzed by the ones of the incident 
waves. 
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Wave vector {2n/a) Transmission z 



Fig. 8 (Left) Transmission spectra of even mode and dispersion diagram of 3D W4 SOI "triangular lattice" CNPW. 
(Right) Field distribution inside W4 SOI CNPW in horizontal (top) and vertical (bottom) planes. Grey levels mark 
electric field amplitude. Black contours correspond to waveguide structure. White dashed lines depict positions of the 
corresponding cuts. Here radius is r = 0.2a, the total nanopillar height is h = 3a, the thickness of Si layer equals to 
a. Dielectric constants of Si and 5*102 were chosen as e = 11.5 and e — 2.1025, respectively. 



The calculated transmission spectrum of a 20 periods long, straight "square-lattice" W4 CNPW is shown 
in the left panels of Fig.|7] There are four modes under the light line as it is shown in the band diagram. In 
Fig. |7]the transmission of the fundamental mode is shown directly together with the dispersion diagram. 
The W4 waveguide displays high transmission efficiency (close to 100%) over a broad spectral range. The 
position of the cut-off frequency is clearly seen in the spectrum. 

In the right panels of figure |7] the band structure and transmission spectra are shown for the W3 
"triangular-lattice" CNPW. A 20 period long, straight CNPW is cut in the F — X direction of the tri- 
angular lattice. A substantial suppression of the transmission is seen in the spectrum, coinciding exactly 
with the position of the mini bandgap in the band structure. Changing the parity of the signal field dis- 
torts the spectrum reflecting the mode symmetries. The even mode displays high transmission efficiency 
(close to 100%) over a broad spectral range. The odd mode has a lower level of transmission and is mostly 
transmitted at higher frequencies. Here, by odd and even modes we understand the corresponding first 
two fundamental modes of a conventional dielectric waveguide. The surprisingly high transmission of 
the even mode above the cutoff frequency, uj « 0.34, can be explained by the resonant behavior of the 
folded radiation mode with negative group velocity [21]. We found similar behavior above cutoff for other 
"triangular-lattice" CNPW structures. 

An example of 3D calculations for an SOI W4 "triangular-lattice" CNPW are presented in Fig. [8] In 
the left panel the dispersion diagram of the structure are presented, while its transmission spectrum for 
the even mode is plotted in the central panel. The dispersion diagram was calculated using 3D supercell 
PWM. In general, the transmission spectrum is very similar to the corresponding spectrum of 2D structure 
(Fig.|7] right). The transmission band is rather broad with 80% transmission efficiency at maximum and a 
sizable stopband at the mini bandgap frequencies. We attribute the moderate level of transmission to the 
impedance mismatch at the conventional waveguide-nanopillar waveguide interface. In the right panel of 
Fig. m the steady-state electric field distribution inside the CNPW is shown for a monochromatic source 
with normalized frequency uj = a/A = 0.3. The field is well confined within the waveguide core in both 
the horizontal and vertical planes, there is no evidence for strong energy leakage into the substrate. 
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Normalized frequency {afk) 

Fig. 9 Normalized energy flux of electromagnetic radiation escaping from the resonator into the terminal (solid 
line) and elsewhere (dashed line) for three nanopillar structures shown in the insets: a Wl with a point defect [24] 
(left); a Wl with Cantor-like longitudinal geometry (center) and a W3 Cantor-like CNPW (right). Here, ds ~ 0.5a, 
cLl = 0.81a, r — 0.15a, A = 0.75a and e = 13.0. The point defect is created by doubling the radius of a periodic 
waveguide with period d^. Arrows mark the resonances discussed in the text. 

4 Aperiodic nanopillar waveguides 

In addition to allowing arbitrary variation of the period and displacement (which is one of the advantages of 
the nanopillar waveguides as opposed to the PCWs), CNPWs allow arbitrary modification of the longitu- 
dinal geometry. A localized change of the properties introduced in one or several nanopillars would create 
a point defect, which functions as a resonator [2,24]. The design of such micro-resonators on the scale of 
a few wavelengths is essential for integrated optics applications. Ideally, such resonators should combine 
the apparently contradictory features of a high Q-factor and of a sufficiently good coupling to a waveguide 
terminal to inject or extract light into or from the resonator Due to the absence of a complete bandgap, 
the breaking of translational symmetry inevitably results in radiation losses of the resonator mode, which 
raises the need for optimizing the Q-factor of the resonator in ID nanopillar waveguides. There have been 
some proposals to decrease the losses based on either mode delocalization [23] or on the effect of multipole 
cancellation [24]. A delocalized mode typically suffers from a decrease of the Q-factor. On the other hand, 
the spatial radiation loss profile of a mode described in Ref. [24] has a nodal line along the waveguide axis, 
which means poor coupling to any components coaxial with the waveguide. 

Other than by means of a point defect, a resonant system can also be created by changing the pe- 
riodic arrangement of nanopillars into a non-periodic one. We show that the use of such aperiodically 
ordered waveguide leads to improved coupling to the coaxial terminal without considerably sacrificing the 
2-factor of the resonant modes. We use fractal Cantor-like NWPs as an example [25]. To construct 
aperiodic NPW, nanopillars of equal radius are aiTanged in a ID chain, where the distances between 
adjacent pillars are given by the Cantor sequence. If we denote S and L for short and long distance 
(ds and d^), respectively, the Cantor sequence is created by the inflation rule L LSL, S SSS and 
unfolds in the following self-similar fashion, which represents a series of middle third Cantor prefractals 
L LSL LSL SSS LSL LSLSSSLSL SSSSSSSSS LSLSSSLSL ^ • • • . 

In order to compare the amount of energy gathered by the coaxial terminal and dissipated elsewhere, we 
excite the system by a dipole source emitting a pulse with a broad spectrum, and use the FDTD method to 
investigate the process of energy loss into the suiToundings. Fig.|9]shows the results. For the point-defect 
structure, the radiation of the resonant mode primarily escapes sideways (Fig. |9] (a)), so despite having 
a high Q-factor (2.1 x 10"*), the coupling between the resonator and other components cannot be made 
efficient. The Cantor structure shows a considerably improved coupling (Fig. |9] (b)) accompanied by a 
drop of the Q-factor down to 2.7 x lO^^. Using a W3 CNPW with Cantor geometry (Fig.|9](c)) raises it 
back to Q 1.1 2.4 x 10'' while still providingjust as good coupling to the coaxial terminal. 

One should notice that the Cantor geometry is only one kind of deterministically aperiodic aiTange- 
ment. Other kinds, e.g., quasi-periodic Fibonacci-like one, can be used leading to a modification of the 
mode structure of NPWs as well as the coupling efficiency of resonant mode into coaxial terminal [25]. 
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Fig. 10 Dispersion diagram for the W2 CNPW section of a directional coupler (left). Directional coupler based on 
two Wl CNPWs (right). Grey levels mark field intensity. Here e = 13.0 and r = 0.15a. 



Engineering the longitudinal geometry of CNPWs appears to be a promising and powerful tool for a further 
degree of freedom in controlling the dispersion properties. 

5 Applications 

Relatively high transmission efficiency and flexibility in dispersion tuning of CNPWs may initiate their use 
as components for efficient and compact nanophotonic devices. Here we discuss two possible applications 
of CNPWs in integrated optics: a coupled nanopiUar waveguide directional coupler [22] and a switchable 
coupled mode laser [31]. 

5 . 1 Directional coupler 

A pair of CNPWs can be used as an effective directional coupler [22]. An example of such a direc- 
tional coupler based on two Wl waveguides is shown in Fig. [10] Analyzing the dispersion diagram of 
the coupling section, namely W2 CNPW (Fig.[TO] left panel), one can see a pronounced difference in the 
propagation constants of the even and odd supermodes in the frequency region around u; = 0.25 — 0.27. 
It is a result of the strong interaction of coupled waveguides, which now are much closer to each other 
than in the case of standard line defect waveguides in a PhC lattice (see, for example, the similar rod 
structure in [26]). In this frequency range the difference between the even and odd supermode propaga- 
tion constant is close to 0.1 • 27r/a, which leads to a crude estimate of minimum coupling length [27]: 
L — TT/\keven — kodd\ = (t^/O.!) • (a/27r) — 5a. Enhanced interaction leads to a shorter coupling length. 
This is illustrated in figure (TO] right panel, where the time averaged squared electric field pattern is shown 
for the normalized frequency ui — 0.26. Guided light hops from the bottom Wl waveguide to the top one 
and back on a distance equal to approximately 5a, which represents well the estimated value. In contrast 
to a directional coupler proposed in [28], the CNPW structure does not require a specially adjusted sepa- 
ration layer between coupled waveguides, thus considerably simplifying a directional coupler design and 
fabrication. 

There are several parameters, which can be used to optimize the directional coupler, e.g., length of the 
coupling region or the number of rows in each of the waveguides. As it has been shown in Sec. |2] lon- 
gitudinal and transverse offsets between the individual waveguides, as well as variation of the dielectric 
constant and radius of the rods substantially modify the dispersion of the compound system, thus affecting 
the coupling efficiency. For example, by shortening the distance between two waveguides one can dra- 
matically increase the propagation constant difference and reduce the coupler size. A similar effect cannot 
be achieved with standard PhC waveguides without any special design tricks involving intermediate walls, 
which increases the complexity of the fabrication procedure. An arbitrary longitudinal offset breaks the 
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symmetry of the device with respect to the symmetry plane between the two Wl waveguides, which may 
further improve the coupling strength similar to the case of an antisymmetric grating coupler [29, 30]. 

5.2 Laser resonators 

The periodicity of the coupled nanopillar waveguides can ensure the distributed feedback within a finite 
waveguide section. This can be seen from the flat tails of the nanopillar waveguide modes near the IBZ 
edge, which correspond to a very low group velocity of CNPW modes. Taking into account that any 
of the CNPW modes can be efficiently excited in the waveguide using an external seeding signal of the 
appropriate spatial profile [14], we have proposed the design of a switchable laser resonator [32] with 
distributed feedback based on the CNPW. The possibility to tune the number of modes, their frequency 
and separation (Sec. |2|i would make such a resonator a promising candidate for a chip-integrated laser 
source. 

The concept of switchable lasing was originally proposed in Ref. [14] and has got further justifications 
in our recent work [32]. In essence, a switchable microlaser comprised of a multimode microresonator, 
where lasing can be switched on demand to any of its eigenmodes by injection seeding [33, 34], i.e. by 
injecting an appropriate pulse before and during the onset of lasing, such that the stimulated emission 
builds up in a designated mode selected by this seeding field rather than from the random noise present in 
the system due to quantum fluctuations and spontaneous emission [32]. 

To provide a basic physical picture of switchable lasing we first consider briefly a simple semi-classical 
laser model in the case of two identical coupled single-mode cavities [32]. In this case there are two 
modes, the symmetric and the antisymmetric one, characterized by spatial field distributions ui,2(r) and 
frequencies wi 2 — ^^o T respectively. Here Auj is the mode detuning from the frequency of the 
single-cavity resonance, luq. For weak mode overlap the spatial intensity profiles of the two modes nearly 
coincide, |Mi(r)|^ « |u2(r)|^. We assume that the cavities contain a laser medium with a homogeneously 
broadened gain line of width Auia > Alo, centered at frequency LOa ~ luq + 5. Here 5 is the detuning of 
the gain profile from the cavity frequency tJo- For this system the semiclassical Maxwell-Bloch equations 
[33, 35], in the rotating-wave and the slowly varying envelope approximation read 

S) 

gR2vC2 (^[allCi~allReix2M2i)] \E,\^ + allC2 \E2\'') E2{t) + F2{t). (9) 

Here all the spatial dependencies of the electric field and atomic polarization were represented in the 
basis of the two cavity modes, such that E{y, t) = £'i(i)Mi(r)e^"^i* + £'2(t)w2(r)e~*'^^*, etc., and the 
atomic polarization was eliminated adiabatically [36,37]. Ej (t) are slowly varying envelopes of two modes 
j = 1, 2. In Eqs. (|9]l the terms linear in Ej{t) describe stimulated emission driving, where the light-matter 
coupling constant is denoted by 5 ~ y/2TrLLjQ(P /h, the pumping rates projected onto the two resonator 
modes byi?j — u*(r)wj(r)i?(r)(ir, and the cavity mode decay rates by Kj. Here d is the dipole moment 
of the atomic transition. The coefficients Cj = Re /3~^, with /3i 2 — Auja/2 + i{S i: Alo), account for the 
different mode-to-gain couplings due to asymmetrical detuning of the atomic transition with respect to the 
resonator frequencies. The terms cubic in Ej (t) describe field saturation above the lasing threshold, where 
7] = (P /2j^^h'^ and the overlap integrals a]^^ = u* {r)uj {r)u'l{r)ui{r)dr are taken over the regions G 
containing the gain medium. Here 7|| is the non-radiative decay rate. The frequency dependence of the 
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Fig. 11 (Left) Cavity phase diagrams for a lasing system governed by Eqs. ^ for FjF{t) ^ Ff{t) in the case of 
symmetric distribution of mode frequencies with respect to gain. The dots denote the stable cavity states and the curves 
represent the phase trajectories for their temporal evolution for different ratios Fl : Fi in the direction of the arrows. 
(Right) Time evolution of the total laser field (shown in grey) and the two cavity modes envelops \Ei \ (solid line ) and 
\E2 1 (dashed line) for the ratio Fl : — 4 : 3. One can see how the first mode wins the competition. Here numerical 
values of the coefficients in Eqs. ([9} were calculated for W2 CNPW resonator with the following parameters, period 
a, d = 1.21a, r — 0.15a and e — 13.0. 



cross-saturation terms is given by A4ij = ^ + (/3*) ^ j- Since |ui(r)|^ « |u2(7')|^ we can further 
assume that a" = a^- = a, Ri = R2 = R and ki — K2 ^ k. 

The inhomogeneous terms Fj (t) originate from the external injection seeding field and from a noise field 
accounting for spontaneous emission [36]. For vanishing functions Fj, Eqs. ^ would take the form of the 
standard two-mode competition equations [33,35], describing bistable lasing [38] and mode hopping in the 
presence of stochastic noise in the system [39]. If both an external seeding field £'^{r,t) and a stochastic 
noise field f "(r, t) are present in the cavity, £{r, t) — £^{r, t) + f "(r, t), the inhomogeneous terms are 
given by, 

Fj{t)^'^^ f dt'e"^^*' f Uj{r)£{r,t')dr = FJF{t)+F^{t). (10) 

T Jt-T Jg 

The time integration in Eq. ( fTOb is the averaging over a time interval larger than 1 / Auj. The function F{t) 
is determined by the temporal dependence of the seeding signal £'*(r, t). The coefficients F^ and F^^{t) 
are determined by the spatial overlap of each mode with the seeding and noise fields, respectively. 

We consider the situation when the seeding prevails over the noise, i.e., FJF(i) ^ -^j^i'^)' before and 
during the onset of lasing. After the onset the Ej become so large that the terms Fj have no effect anymore. 
In this situation the evolution of the resonator will be determined by the ratio of F^ and F|. In Fig.fTTIdeft) 
the phase trajectories of the temporal resonator state evolution in the (|£'i | , [£'2 1 ) space is presented for 
different values of F^ and F|. As seen in Fig. [TT](left), the lasing state first reaches overall intensity 
saturation (|£'i| + \E2\ = E^) and then drifts towards one of the stable fixed points corresponding to 
single-mode lasing (either = E'^ or |i?2| = E^). The drift happens on a longer time scale than the 
initial overall intensity growth, and the intermode beats decay fast after the lasing onset (Fig. [TT] right). 
The drift occurs towards the mode whose spatial and temporal overlap with the seeding signal is larger, 
demonstrating a switchable lasing behavior. It is important to note that even in the case of asymmetric 
detuning of the cavity modes with respect to the gain frequency (5 7^ 0), single-mode lasing is achieved 
into the mode whose spatial overlap with the seeding field, F^, is largest, i.e., if one of the following 
conditions, Ff ^ or Ff ^ F|, is satisfied. 

To demonstrate the predictions of this simple theory we have modeled the lasing action in four-row 
CNPW structure with a realistic injection seeding (Fig.fT2li using the FDTD method [31]. The externally 
pumped laser-active medium is placed in the central 7 pillars of all four rows. This is done to maximize 
coupling between the active medium and the main localization region of the lasing modes. The population 
dynamics of an active medium is described at each space point by the rate equations of a four-level laser 
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Fig. 12 The amplitude spectra (top) and laser filed distribution (bottom) for the periodic injection- seeded four-row 
CNPWs. The lines labeled "Mode-l" to "Mode-4" correspond to the seeding signals (a)-(d) shown in bottom panel. 
The shaded areas represent the laser amplification line, with its central frequency uja = 0.3225. The pumping rate 



equals Wj, = 1.0 x 10' 
excited in the terminal. 



The panels (a)-(d) correspond to different seeding signals, shown schematically as 



with an external pumping rate Wp. To achieve population inversion we have chosen the following values 
for the non-radiative transition times, T32 ~ tiq <C T21, with t^i = rio = 1 x 10^ ^'^ s, T21 = 3 x 10~^° s, 
and the total level population is iVtotai = 10^'' per unit cell [40]. The Maxwell equations are solved using 
FDTD scheme supplemented by the usual equation of motion for the polarization density in the medium 
and by the laser rate equations [40^3]. All calculations were done for TM polarization. The seeding 
signal is excited by four emitters (linear groups of dipoles) engineered on the regular dielectric waveguide 
attached to the CNPW structure (see Fig.fTSI). Each of the emitters generates a single short Gaussian pulse 
with carrier frequency LOa and with half-width duration at = lO'^dt. The relative phase of the fields in 
these pulses is chosen or tt. Technically, the seeding dipoles are realized as point like oscillating current 
sources in the Maxwell equations [31]. Similarly, the spontaneous emission [42,44,45] can be modeled 
as an ensemble of point current sources, randomly placed in space, with temporally (5-correlated Langevin 
noise [45]. The computational domain of size 7a x 22o was discretized with a mesh point spacing of a/16. 
The time step is related to the spatial mesh to assure stability and was chosen dt = Qx 10^^^ s. To simulate 
an open system, perfectly matched layer (PML) boundary conditions [20] were used. 

In Fig.[T2](top) the lasing spectra in the steady state long after the seeding signal has decayed is shown. 
The broad shaded areas depict the laser line of width AcJq centered at uJa, which is shifted slightly towards 
lower frequencies. As a rule the Q-factor is larger for modes with the higher frequency. The shifted 
laser line compensates this Q-factor difference, so that any of the four CNPW modes can be selected by 
the appropriate seeding signal with the same symmetry. In figure [12] (bottom) the spatial electric field 
distribution in the four-row CNPW laser resonator is shown at an instant of time long after the seeding 
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signal has decayed and after the steady state has been reached. The symmetry of the selected lasing modes 
corresponds to that of the seeding signal (Fig.[T2]l. 

The proposed concept of switchable lasing is not limited to the periodic CNPW structures, but is ex- 
pected to work in any resonator featuring bi- or multistability. Any coupled cavity based system would be 
a good candidate for the effects predicted. For example aperiodic CNPW based resonator also show the 
switchable lasing behavior for resonant modes discussed in Sec.|4][31]. 

6 Conclusion 

We have shown that a novel type of coupled nanopillar waveguides, comprised of several periodic or 
aperiodic rows of dielectric rods, may have potential applications in compact photonics. The strong cou- 
pling regime can be utilized in ultrashort directional couplers or laser cavities, which might possess an 
additional functionality and flexibility when different longitudinal and transverse offsets among individual 
waveguides are employed. The factors of major influence upon the mode dispersion have been analyzed. 
Transmission spectra for 2D and 3D systems prove the possible single mode excitation by imposing spe- 
cific symmetry conditions onto a field source and high transmission characteristics of coupled nanopillar 
waveguides. 
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